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ABSTRACT 

Using a set of high-resolution spectra, we studied the physical and orbital 
properties of the O-type binary CPD — 41°7733, located in the core of NGC6231. 
We report the unambiguous detection of the secondary spectral signature and we 
derive the first SB2 orbital solution of the system. The period is 5.6815±0.0015 d 
and the orbit has no significant eccentricity. CPD — 41°7733 probably consists of 
stars of spectral types 08.5 and B3. As for other objects in the cluster, we ob- 
serve discrepant luminosity classifications while using spectroscopic or brightness 
criteria. Still, the present analysis suggests that both components display phys- 
ical parameters close to those of typical 08.5 and B3 dwarfs. We also analyze 
the X-ray light curves and spectra obtained during six 30 ks XMM-Newton point- 
ings spread over the 5.7 d period. We find no significant variability between the 
different pointings, nor within the individual observations. The CPD — 41°7733 
X-ray spectrum is well reproduced by a three-temperature thermal mekal model 
with temperatures of 0.3, 0.8 and 2.4 keV. No X-ray overluminosity, resulting e.g. 
from a possible wind interaction, is observed. The emission of CPD — 41°7733 is 
thus very representative of typical O-type star X-ray emission. 

Subject headings: stars: individual: CPD — 41°7733 - stars: binaries: close - 
stars: binaries: spectroscopic - stars: early-type - X-rays: individuals: CPD — 41°7733 
- X-rays: star 
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1. Introduction 

Early-type stars of spectral type O, and their evolved descendants, the Wolf-Rayet 
stars, are among the hottest and brightest objects in the Galaxy. Through their large 
radiative and kinetic energy input, they have a strong impact on their surroundings and 
their study is relevant for a number of galactic and extra-galactic issues. However, despite 
their importance, a proper mapping of their physical parameters is still lacking. This mainly 
results from the limited amount of observational constraints available so far. Indeed, only 
about 30 d ouble- lined spectroscopic binaries (SB2) with O-type components have definite 



orbits (e.g. iGiesI 120031 ) . In this context, the present paper is devoted to a rather poorly 



known object, CPD — 41°7733, which lies in the core of the young open cluster NGC 6231. 



Sune et al. 


1998; 09 III, 


Walborn 


1973) 



(119441 ) because of the large difference between the star velocity and the cluster mean velocity. 
Since then, three SB1 orbital solutions have been published, but yielded discrepant orbital 
parameters (see §3.31) . Clearly, a new study, based on high quality data, was needed to 
search for the companion signature and to definitely upgrade the orbital (§3]) and physical 
O parameters of this system. In we also investigate the CPD — 41°7733 X-ray properties 
using a recent 180 ksec XMM-Newton observation of the NGC 6231 cluster. Finally, we give, 
in §nj a summary of the main results of this work. 



2. Optical and X-ray observations 

The present work is based on over 30 high-resolution high S/N spectra of CPD — 41°7733 
obtained at the European Southern Observatory (ESO, La Silla, Chile) and spanning a 
time interval of six years. Table [T] provides the Julian dates at mid-exposure time, the 
instrumental setups and the mean radial velocity (RV) at each epoch. Prior to averaging, 
the individual line RVs were referred to a zero-systemic velocity frame using the systemic 
velocities quoted in Table [3] (the considered lines are listed in Table [3]). For some epochs, 



1 The SI MBAD astronomical data base reports the ISeggewisd number of CPD — 41°7733 to be 297 
According to iRaboud et al.1 (|1997h , this number belongs t o an extension of thelSeggewisd (119681) catalogue - 



1992). However, 



Raboud et al 



that o nly listed 295 stars - included in WEB DA data base lMermilliodl |l988l 
([19971 ) quoted a magnitude of 12.5 for the star Se 297. This does obviously not correspond to CPD — 41°7733. 
From the positions on the CC D, the magnitudes a nd the spectral t ypes o f the objects listed in their catalogue, 



we found that the 
Se323. 



Seggewisd number reported bv lRaboud et al.l (|l997j ) for CPD — 41°7733 should rather be 
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no reliable measurement could be obtained for the secondary, either because of the low S/N 
ratio or because of the important blend between primary and secondary lines. In those cases, 
the corresponding col umn is left empty. Instrumental setups and reduction techniques have 
been described in e.g. ISana et al.l (l2006bl ) and will not be repeated here. 



In parallel to the optical monitoring, six 30 ks X-ray exposures of the NGC 6231 clus- 
ter were performed by XMM-Newton in September 2001. The campaign was described in 
Sana et al.l (j2006al ) and we restrain the description here to additional elements specific to 
CPD — 41°7733. For each of the six pointings and for the three epic instruments, we ex- 
tracted broad band X-ray light curves and spectra using the SAS task evselect. The circular 
extraction region was centered on the source position and, due to the presence of three 
bright X-ray neighbors (HD 152249, HD 326329 and NGC 6231 SBL 324), was limited to a 
radius of 1 1'/5 (Fig. [T]). The different backgrounds were estimated from th e very few source- 
free regions in the field of view (for more details, see ISana et al.ll2006d ). Finally, we also 
extracted the merged spectra for each instrument, thus combining the six observations of 
CPD — 41°7733. For this purpose, we built the corresponding arf files using the SAS task 
arfgen, and we adopted the rmf files provided by the SOC and adapted to the particular 
position of CPD — 41°7733 on the detectors. Julian Date (JD) at mid-exposure, effective ex- 
posure times and background-corrected count rates as obtained for the different instruments 
are quoted in Table El 



3. CPD -41°7733 orbital solution 



3.1. Optical spectrum 



The spectrum of CPD — 41°7733 (Fig. [2]) is clearly dominated by the H I Balmer, He I 
and He II absorption lines. Numerous metallic lines (C, N, O, Si and Mg ions) can also 
be identified. The Balmer lines and all the He I lines present a clear SB2 signature with 
the primary lines being several times stronger than the associated secondary lines (see e.g. 
Fig. [3]). In the very high S/N ratio spectra, the secondary signature might be seen in a 
few metallic lines (e.g. C IIA4267, Si III AA4552-68-74) but this remains at the very limit 
of detection. Two emission lines, S IVA4486 and C IIIA5696, are further detected in the 
CPD — 41°7733 spectrum; both are associated with the primary component. 

Line positions and equivalent widths (W\) were measured by adjusting G a ussian curves 
to the studied pro files. Effective w avelengths for O-stars from IConti et al.l (119771 ) below 
4800 A and from lUnderhilll ( 19941 ) above were adopted to compute the radial velocities 
(RVs). For the metallic lines that are not listed in these latter works, we used the rest 
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wavelengths from iMoord (119591 ) 



3.2. Orbital solution 



To search for the orbital period P of the s ystem, we applied both th e method of 



2001 



Lafler fc Kinmanl (119651 ) and the Fourier analysis of iHeck et al.l (119851 . see also lGosset et al. 



for comments) on the RV sets associated with the He IA4471, He IIA4686, O III A5592 
and He IA5876 lines. The obtained periods are all in excellent agreement with an average 
value of 5.68156 d. Estimating the period uncertainty on the basis of one tenth of the peri- 
odogram peak width, we found op = 1.5 to 1.8 10 _3 d according to the data set considered. 

We then computed orbital solutions using individual RV sets associated with the differ- 
ent absorption lines listed in Table [3j For SB2 solutions, the RV equations were transformed 
in order to sym metrically propagat e the errors on the parameters associated with each of 
the components (ISana et al.l l2006bl ) . We adopted a period value and a relative primary to 
secondary weight ratio (s y /s x ) that yield the lowest \ 2 - We tested both circular and ec- 
centric orbits but found that the latter did not improve the quality of the fit. Table [3] lists 
the orbital solutions obtained for the data sets associated with different lines. They are all 
in excellent agreement with each other. As a final step, we computed the averaged RVs of 
all the primary lines and of the SB2 He I lines. For this, we shifted the individual RVs 
to a common velocity frame, taking into account the different systemic velocities deduced 
from the individual orbital solutions quoted in Table [31 The orbital solutions obtained using 
either the averaged primary line RVs only or the averaged He I SB2 RVs are given in Table H] 
(respectively labeled 'Prim.' or 'He I lines'). The quoted period uncertainties correspond to 
the 1-cr confidence interval^. We emphasize however that the uncertainty on the periodicity 
of the phenomenon (which is independent of any model consideration) is rather related to 
the peak width in the periodogram as stated at the beginning of this section. The other 
uncertainties (always quoted as one standard deviation) were computed by error propagation 
from the least-square fit. Finally, we adopted the apparent systemic velocities of both com- 
ponents as the weighted means of the systemic velocities in the different He I-line solutions 
of Table [3J 7! = -21.5 ± 0.2 km s" 1 and ~ 2 = -23.1 ± 0.2 km s" 1 . 



2 The confidence intervals on the period value were computed by varying the period (all other parameters 
being kept constant at their best-fit value) until the minimum \ 2 value of the fit is increased by a value 
A% 2 . The latter value was chosen to correspond to the 68.27% confidence level and depends on the number 
of degrees of freedom of the fit. 
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3.3. Towards a global solution including published and new data 



As mentioned earlier, three SB1 orbital solutions have previously bee n published. Based 



Perry et all (11990 



HCB74I hereafter) 



on a set o f 16 observations spread over 8.1 days plus four observations from 
PHYB90l hereafter) obtained about 800 days earlier, IffiU et al.l (11974 
derived the first SB1 orbital solution with a period P = 5.64±0.01 d, a small eccentricity e = 
0.04 ±0.03 and an amplitude for the primary radial velocity cur ve of = 89.0 ±2.7 km s -1 . 



With two more observati ons and including r adial velocities from lStruvd (I1944T) andlHill et al. 
(11974 ) (but n ot those of iPerry et al.l Il990l . only published much latter) , iLevato fc Morrell 



(119831 . ILM83I hereafter) derived a new orbital solution. They found a larger period (P 
5.74973 ± 310 _5 d) and eccentricity (e = 0.14 ± 0.015) while the amplitude of the radial 
velocity curve remained mostly unchanged (K\ = 84 ± 2kms~ 1 ) . Finally, CPD — 4 1 °7733 



was o bserved once by IUE. Including this additional measurement, IStickland &: Lloydl (12001 
SLOll hereafter) published a slightly modified orbit with P = 5.749809 ± 3.010 _5 d, e = 
0.0 45 ± 0.021 and Ki = 83.9 ± 2.1 km s -1 . They stated however that they could not include 
the |PHYB90l data in their fit. 



Though the present data set yields orbital parameters of the same order of magnitude 
than the ones obtained from previous works, fundamental contradictions exist between the 
different studies. Our final period P = 5.681504 d is significantl y differ ent from the previous 
determinations and we report a circular orbit while, for example, ILM83I quoted an eccentricity 
up to 0.14±0.015. Our derived semi-amplitude of the primary RV curve K\ is systematically 
higher tha n the previously published values, with differences up to lOkms -1 compared to 
LM83I and lSLOll . It was thus crucial to check the previous orbital solutions using the different 
RV sets adop ted by these authors. Doing this, we noted that the second m easurement listed 
bv |PHYB90l at RV = -109 km s" 1 , is clearly not fitted by the 1HCB74J solution though it 
was included in the adjustment. Actually its location in the p hase diagram is shifted by 
about 0.2 in phase compared to the figure displayed by IHCB74 Thi s prob ably results from 
a typogra phic error in the published Julian date, which explains why lSLOll could not include 
PHYB90 observations in their combined solution. 



Though we have been able to reproduce the results of IHCB74I and jSLOll using their 
respective data sets, we could not reproduce the results proposed by lLM83l . Both the period 
value and, more particularly, the eccentricity derived from their data differ significantly. In 
all our computations, whatever the data set used, the maximum eccentricity that we found 
was e = 0.07, with a null eccentricit y well within a 2-cr interval. W e thus defin itively rule 
out the large eccentricity quoted by ILM83I . Furthermore, both the ILM83I and ISLOll data 
sets actually carried a strong ambiguity on the period, with more than ten strong aliases 
between 5.6 and 5.8 days. Clearly, their period determination was far less constrained than 
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what their error bars suggested. 



Finally, we combined all the available primary RV measurements (except the IPHYB90 
point at RV = — 109 km s -1 ) to derive a combined orbital solution. For the RVs from the 
present work, we used the averaged RVs reported in Tabled] to which we added the adopted 
primary systemic velocity 71 = — 21.5 km s -1 . A period search using both the Lafler & 
Kinman and the Fourier analysis techniques (see §3.21) yield values close to 5.6815 d with an 
uncertainty of o~p = 2.4 10~ 4 d. Again no significant improvement of the quality of the fit 
is obtained assuming an eccentric orbit. Fig. [5] shows the values of the fit r.m.s. for periods 
ranging from 5.6 to 5.8 d. Clearly, the present analysis allows us to solve the ambiguity 
about the period that plagued earlier investigations. The best-fit orbital parameters are 
listed in Table H] (labeled 'Lit.') and the RV curve is displayed in Fig. [6J The corresponding 
heliocentric systemic velocity is —21.7 ±0.1 kms -1 . 



4. CPD -41° 7733 physical parameters 
4.1. Spectral types and luminosity classes 



4.1.1. Primary component 



The spectral signature of the primary component is easily discerni ble in the spectrum of 
CPD —4 1 7733 . Adopt ing the classification criteria from lContil (Il973bl ) as adapted to late O- 
stars by iMathyj (1988), we obtained a mean logVy'd^^) = 0.22 ±0.05 which corresponds 
to a spectral type 08.5, with spect ral type 08 within 1 a. To determine the luminosity 
class, we adopted the criterion from lConti k, Alschulerl (Il97ll ). We obtain log W(^^-) = 
0.27±0.04, which leads to a giant lumino sity clas s . We also measured log W" = log(W / A438s)± 
log(WA4686) = 5.23 ±0.04. According to iMathysl (119881 ). this rules out a supergiant class and 



points towards a giant classification, unless the ratio l\ = -j^ < 0.67. The optical brightness 
of both components will be discussed in §4.21 



4-1-2. Secondary component 



The main spectral signatures of the secondary star in CPD— 41°7733 are the Balmer 
and He I lines. The absence of the He II lines, as well as of the O III A5592 line, at the 
positions predicted by the orbital solution definitively excludes an O spectral type and, 
at our detection threshold, m ost probably indicates a spectral subtype later than B0. 7 
(IWalborn fc Fitzpatricklll990l ). To refine our classification, we searched for the presence of 
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secondary metallic lines. This led us to consider low intensity lines, with equivalent widths 
(relative to the composite spectrum continuum) down to ~0.0lA. We therefore focused on 
the six FEROS spectra obtained at the ESO/MPG 2.2m telescope, that exhibit S/N ratios 
above 200. Mean primary and secondary equivalent widths as measured on these spectra 
are reported in Table [51 Some of the secondary lines are clearly deblended from neighboring 
lines, and thus securely measured, in only one or two of these six spectra. These are marked 
with a colon (:) and should be considered as a rough indication of the line strength only. 

The absence of He IIA4686, C IIIA4650, N IIIA4097 and Si IVAA4089, 4116 lines in 
the secondary spectrum, combined with the presence of a comparatively strong C IIA4267 
line, point towards a spectral type later than B2 while the lack of Fe II lines corresponds 
to a spectral type earlier than B5. The usual luminosity criterion in this spectral range is 
Si IIIA4552/He IA4388. From our measurements, we obtained a ratio of about one to ten; 
this result definitively excludes a supergiant luminosity class, which is in agreement with the 
low intensity of the secondary O II spectrum. Similarly, this ratio seems to indicate a main 
sequence luminosity class rather than a giant class. 



4.2. Optical brightness ratio 



We estimated the optical brightness ratio based on the dilution of the primary lines. 
We compared their mean W\ with typ ical (averaged) W\ of 08.5 III stars (IConti fe Alschuler 
197ll : Icbntll973al : lMathvslll988l . ll989h . Based on the He I AA4026, 4388, 4471 and He II AA4542 
4686 lines, we obtained an averaged brightness ratio corresponding to U = 0-8 7 ± 0.10. We 



then compared the secondary intrinsic W\ with typical W\ for B stars (jDidelon!ll982l ) using 
a dilution factor l 2 = 0.13. The intensity of the detected lines in the secondary spectra is in 
rough agreement with a B3 III-V classification but suggests a slightly larger luminosity ratio 
in favor of the primary component. Adopting l\ = 0.87 ±0.10 yields log W"^ ~ 5.29 ±0.06, 
which corresponds to a giant luminosity class for the primary though we note that the 
He IIA4686 line is stronger than in typical 08.5 giants. Should we assume that the primary 
is a main sequence star, the comparison of the observed W\ with typical W\ for 08.5 V stars 
yields h = 0.74 ± 0.08. This yields \ogW^ im « 5.36 ± 0.06 thus in contradiction with the V 
luminosity class hypothesis. Clearly the spectroscopic classification criteria points towards 
the primary being a giant. 
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4.3. Luminosities and stellar radii 



Typical absolute visual ma gnitudes for 08.5 stars are My = —5.2 and My = —4.5 /— 4.4 
for class III and V respectively ([Humphreys fc McElroylll984l ; iHowarth &: Prinjalll989l ). The 
absolute visual magnitude s of B3 stars are My = —3.0 and —1.6 for gian ts and main- 
sequence stars respectively ( jSchmidt-Kalerlll982l ; iHumphreys &: McElroylll984j ). To compute 
the visual magnitudes of the CPD — 41°7733 components, we adopted the cluster distance 
modulus DM = If .07±0.04 as obta ined from the average of the different photometric results 
since the 1990's (feana et al.l l2006ah and we used V = 7.875, R = 3.3 and (B -V) = 0.158 
from ISung et al.l (119981 ) . Using l\ = 0.87 ± 0.10 as derived in the previous paragraph in 
the case of a giant primary and assuming that the secondary is also a giant, we obtained 
Mv = —4.67 ± 0.15 for the total magnitude of the system, and My,i = —4.52 ± 0.20 and 
My t 2 = —2.46 ± 1.12 for the primary and secondary respectively. Under these hypotheses, 
the primary is clearly fainter than typical 8.5 giants. Adopting th e effec tive temperature 
calibration and the bolometric correction of IHumphreys &: McElroyl (119841 ). we constrained 
the radii to values of Ri = 8.8 ± 1.0 R and R 2 = 6.3 ± 3.9 R Q . Again, the stars are too 
small compared to typical giants (R° 8 - 5 111 « 13 R Q and R B3 111 9 R ). 

Therefore, if we assume that the CPD —41° 7733 components are indeed main sequence 
objects, the same reasoning with l\ = 0.74 ± 0.08 yields M v = —4.64 ± 0.13, M Vi i = 
-4.31 ± 0.17 and M v , 2 = -3.18 ± 0.37, log(^) = 4.98 ± 0.08, log(^) = 3.89 ± 0.17 
and therefore R\ = 8.4 ± 1.0 R Q and R 2 = 8.4 ± 1.9 R . Corresponding typical radii are, 
in this case, of 9.0 R and 4.2 R for 08.5 and B3 main sequence stars. We now observe a 
better agreement for the primary, but the secondary parameters are at odds with those of 
typical B3 dwarfs. 

The best agreement between the deduced and typical parameters is obtained assuming 
that CPD — 41°7733 harbors an 08.5 V and a B3 III component. In the latter case, we derived 
M v ,i = -4.31 ± 0.17, M v ,2 = -3.18 ± 0.37, R x = 8.4 ± 1.0 R and R 2 = 8.7 ± 1.7 R . This 
option is however rather unlikely. We finally note that a larger value for l\ would significantly 
decrease the obtained secondary radius without affecting much R 1 . Such a ratio largely in 
favor of the primary (li > 0.9) is indeed suggested by the typical 08.5 and B3 luminosities. 

In summary, though the spectral criteria rather indicate a giant luminosity class for the 
primary, the estimated values for its magnitude and radius are more consistent with the pri- 
mary being a main sequence star. This situation is reminiscent of the case of CPD— 41°7742, 
an eclipsing SB2 binary in NGC6 231 for which the spectroscopic criteria clearly indicated 
a giant class. In ISana et al.l (120031 ) . we inferred, from the system luminosity, much smaller 
radii than expected for giant stars. This was confirmed by the analysis of the system light 
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curve (jSana et al.l 120051 ) and the resulting parameters are in good agreement with the typ- 
ical values expected for dwarfs of the corresponding spectral types. The locations of the 
CPD — 41°7733 components in the H-R diagram (Fig. |7j) indicate a primary evolutionary age 
of 3 to 4 Myr, in good agreement with previous age estimates for the cluster. This suggests 
that the primary is only just leaving the main sequence and could explain why the star 
displays properties somewhat intermediate between those of a dwarf and a giant. In this 
respect, the estimate of the brightness ratio based on the dilution of the primary lines might 
have been biased. It is thus probable that both components display stellar radii typical of 
main sequence objects. 



4.4. Masses and orbital inclination 



Typical masses for 08.5 V (resp. Ill) stars are about 19 (resp. 24) M (IMartins et al. 



20051 ). Comparing this with the minimal values obtained in Table H] indicates that the orbital 
inclination of the system is probably around 74° (resp. 63°). This corresponds to a secondary 
mass around 7 to 9 M fll which is slightly lower than typical values of B3 III stars (see e.g. 
Schmidt-Kalerlll982l ). but is consistent with B3 main sequence stars. Adopting the different 
values for the radii derived in the previous section, we found that CPD — 41°7733 should not 
display eclipses un less i > 65 — 70° . Clearly th is is a limiting case. Referring to a paper 
in preparation by iPerry et al.l . iHill et al.l (119741 ) announced that the star was displaying 
magnitude variations larger than O^ IS. This re s ult wa s however not confirmed once the 
mentioned pap er finally appear ed as IPerry et al.l (Il990l ) . From our photometric campaign 
of the cluster (jSana et al.l 120051 ) . the r.m.s. of the data set associated to CPD — 41°7733 is 
about CP012, slightly larger than the expected noise. The peak-to-peak variations over the 
one month time span of our campaign is 0.05 mag at maximum. A power spectral analysis 
indicates a slightly dominant peak at v ph 0.084 d _1 in both filters, thus tentatively suggesting 
a period close to 12 d although this can not be the signature of eclipses. 

As a last point, we estimated the projected rotational velocity vsini of the primary star 
by comparing the full widths at half maximum of several lines with those given by model 
spectra computed with an effective temperature and a gravity value corresponding to the 
above estimates for the primary component^]. Based on the lines quoted in Table [3] (with 
the exception of He I AA4471 and 7065), we derived vsini pz 83±8kms _1 . Given a probable 



3 The model spectra were computed, for a grid of -usini ranging from 50 to 300 km s -1 with a step of 
lOkms -1 , using the tlusty and SYNSPEC codes ( Lanz fc Hubenv 20031 ) that use line blanketed, NLTE, 
plane-parallel, hydrostatic atmospheres. 
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inclination of 65-75°, this value is slightly larger, still compatible within the error bars, with 
a synchronous rotation rate (corresponding, for Ri = 8.4 — 8.8 R , to v ~ 75 — 78kms~ 1 ). 
No significant difference was found assuming either a dwarf or a giant luminosity class for 
the primary. 



5. X-ray analysis 



We first extracted broad band X-ray light curves of CPD — 41°7733 using the average 
count rates during each observation. The different energy bands considered are the total 
band [0.5-10.0 keV], a soft (S x ) band [0.5-1.0 keV], an intermediate (M x ) band [1.0-2.5 keV] 
and a hard (Hx) band [2.5-10.0 keV]. We found no consistent variations between the three 
instruments. A \ 2 test performed using the count rates in the different energy bands did not 
allow us to reject the null hypothesis of a constant count rate during all six pointings, except 
for the pn instrument that displays a deviating point in the Mx band at ~ 0.18. Using 
the individual event lists obtained during each of the six pointings, we also searched for 
short term variability. We performed a tim e series analysis s imilar to the one carried out for 
HD 152248, the central target of the field (jSana et al.ll2004j ) . Again no significant variation 
could be consistently detected for any of the three instruments. This strongly suggests that, 
at our detection threshold, the X-ray emission from CPD — 41°7733 is mostly constant. 

To constrain the physical properties of the emitting plasma, we adjuste d a series of 
optically thin thermal plasma mekal models (IMewe et al.lll985t lKaastralll992l ) to the spec- 
tra obtained during each XMM-Newton observation. The epic M PS and epic pn spectra 
were adjusted simultaneously using xspec v.11.2.0 (lArnaudlll996l ). We adopted an equiv- 
alent interstellar column of neutral hy drogen of -/Vh.tsm = 0.26 x 10 22 cm -2 , corresponding 
to E(B — V) — 0.447 as obtained from ISung et al.l (119981 ). A single temperature model was 



insufficient to adequately describe the observed spectra. We thus adopted two-temperature 
(2-T) models allowing for possible local absorption for both mekal components. The absorp- 
tion column associated with the lower temperature component tends to be systematically 
close to zero. Lower residuals and more stable solutions are obtained b y fixing this colum n to 
zero. This situa tion is reminiscen t of what was observed for HD 152248 (jSana et al.ll2004l ) and 
CPD-41°7742 JSana et al.ll2005h . two other early-type binaries in the core of the NGC6231 
cluster. Therefore, the actually fitted model is (wabsisM * (mekal! + wabs 2 * mekal 2 )), 
in which the term wabsigM was fixed to the interstellar value (iV*H,isM = 0.26 x 10 22 cm~ 2 ). 
Table [6] provides the best-fit parameters (iVfi, the absorbing column; kT, the plasma com- 
ponent temperature; norm, the normalization factor) and the limits of the 90% confidence 
intervals. Corresponding X-ray fluxes in the total energy band (fx) as well as in the soft 
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(/x,s)> intermediate (/x,m) an d hard (fx,n) bands are also provided. Finally, for each of 
the epic instruments, we extracted the combined spectrum from the merging of the six 
X-ray observations. We also used a 2-T mekal model to fit the obtained spectra and the 
best parameters are given in the last line of Table [6] The spectra and the best fit model 
are presented in Fig. [BJ The ISM absorption corrected fluxes (f un ") computed in the dif- 
ferent energy bands considered here as well as the observed X-ray luminosity in the total 
0.5-10.0 keV band are given in Table [71 

It is clear from Fig. [8] that an additional, higher energy component is present in the 
spectrum of CPD — 41°7733. Adjusting a three-temperature (3-T) mekal model reveals that 
the third component has a best-fit temperature close to 2.4 keV. We also adjusted a 2-T 
mekal + power-law (PL) model. To avoid the fitting procedure of the PL component to be 
biased by small discrepancies at lower energy, we first adjusted the 2-T model in the 0.5-4.0 
keV range. We then held these parameters fixed to their best fit values while extending the 
energy domain up to 10.0 keV. Again, both the 3-T and the 2-T+PL models fit the observed 
spectra adequately. The best-fit photon index of the PL component is about T »s 2.9, but it 
is very poorly constrained. A PL component with T = 1.5 fits the spectra almost equivalently 
well. In NGC 6231, such a higher energy component (with a typical temperature of a few 
keV) is seen in four objects (HD 152248, HD 326329, object [HD 152314]HD 152314 and 
CPD — 41°7733), w hich corresponds to about one fourth of the O-type stars in the cluster 



(jSana et al.ll2006d ). Three of them are binaries but only one, HD 152248, displays a wind- 
wind collision. The X-ray emission from other systems similar to CPD — 41°7733 (such 
as HD 152219) does not display such an additional higher energy component. Differences 
between the systems with and without such a component remain unclear. 

Finally, we used the constraints on the physical parameters deduced in the previous 
sections to get more insight into the winds of the CPD — 41°7733 components. We assumed 
an orbital inclinat ion of 70°. Mass-loss rates and terminal velocities were estimated following 



Vink et al.l (120011 ). As expected, the primary wind is overwhelmingly dominant and no ram 



pressure equilibrium is possible on the system axis. This suggests that t he pr i mary wind 



will crash into the secondary star surface. However, using the formalism of lUsovl (119921 ). one 
estimates that the amount of X-ray emission that could be produced by such an interaction 
is about 10 29 ergs -1 , thus two orders of magnitude smaller than the intrinsic contribution 
of the two stellar components. As the CPD — 41°7733 orbit is circular, we further do not 
expect to observe an intrinsic modulation of the emitted flux resulting from a variation of 
the shock strength, which is thus in agreement with our observations. The X-ray emission of 
CPD —41° 7733 is thus very representative of the X-ray emission of normal O-type stars and, 
indeed, it is very well fitted by the newly ob tained canonical re lation log(-^-) = — 6.91±0.15, 



derived for the O-type stars in NGC 6231 (jSana et al.ll2006d ). 
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6. Conclusions 

We present the results of a high-resolution spectroscopic campaign on the O-type binary 
CPD — 41°7733. We report the first detection of the secondary spectral signature and we 
derive the very first SB2 orbital solution for the system. The orbital period is close to 5.681 d 
and the orbit is most probably circular. The new orbital elements are significantly different 
from those obtained in previous works. Using the same data sets as previous authors, we find 
that their period values were poorly constrained because of a strong aliasing that resulted 
from the spread of a small number of observations over a large time span. We also combined 
all the observations available from the literature in a joint orbital solution. The resulting 
orbital elements are in perfect agreement with those found using our data set alone. 

CPD — 41°7733 most probably consists of an 08.5 star plus a B3 companion. Though 
different criteria yield discrepant luminosity classifications, both stars probably display phys- 
ical parameters close to those of typical 08.5 and B3 dwarfs. Though the orbital inclination 
is rather large (around 65-75°), our photometric campaign on NGC6231 did not reveal 
any eclipse. Finally, we estimated the projected rotational velocity of the primary to be 
around 85kms _1 , slightly larger than, but still compatible within the error bars, with the 
synchronous rotation rate. 

We also analyze the X-ray light curves and spectra of CPD — 41°7733 obtained during 
the six XMM-Newton pointings towards the cluster. We find no significant variability between 
the different pointings, nor within the individual observations. The X-ray spectrum is well 
reproduced by a two-temperature thermal mekal model with k7\ as 0.3 keV and kT 2 R3 
0.8 keV. The merged spectrum built from the combination of the six pointings shows an 
additional higher energy component. The latter can be described either by a mekal model 
with an energy close to 2.4 keV or by a power-law component with a photon index T close 
to 3, though the latter value is very loosely constrained. No X-ray overluminosity resulting 
from a possible wind interaction phenomenon is observed and, indeed, none is expected 
at our detection level. As a consequence, the emission of CPD — 41°7733 should be very 
representative of typical O-type star X-ray emission. 

It is a pleasure to thank the referee, Dr. D. Gies, for detailed and helpful suggestions, 
as well as N. Linder for running the synspec and tlusty codes. The authors are greatly 
indebted towards the FNRS, Belgium. This work was supported by the PRODEX XMM 
and Integral Projects, contracts P4/05 and P5/36 'Pole d'Attraction Interuniversitaire' (Bel- 
gium) . 

Facilities: ESO/Max Plank:2.2m (FEROS), ESO:1.52m (FEROS), ESO:CAT (CES), 
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Table 1: Journal of the spectroscopic observations 

HJD H eI RV ljX - 7 1,A RV 2 ,X - 72,A 



-2 450 000 




(kms 1 ) 


(kms 1 ) 


996. 


684 a 


0.259 


96.8 


-261.6 


997. 


638 a 


0.427 


40.3 


-107.4 


998. 


627 a 


0.601 


-55.1 


149.3 


999. 


634 a 


0.778 


-89.8 




1000. 


607 a 


0.950 


-28.7 




1299. 


827 b 


0.615 


-62.3 


156.6 


1300. 


818 b 


0.790 


-90.3 


238.2 


1301. 


825 b 


0.967 


-14.6 




1302. 


827 b 


0.143 


74.1 


-190.8 


1304. 


822 b 


0.494 


1.0 


3.5 


1327. 


843 b 


0.546 


-24.8 


83.8 


1668. 


865 b 


0.570 


-39.2 


108.5 


1669. 


867 b 


0.746 


-93.6 


246.8 


1670. 


858 b 


0.920 


-45.4 


120.0 


1671. 


862 b 


0.097 


53.7 


-137.6 


1672. 


851 b 


0.271 


92.6 


-247.0 


2037. 


839 b 


0.513 


-4.6 


-1.3 


2039. 


826 b 


0.862 


-71.9 


188.3 


2040. 


837 b 


0.040 


22.4 


-55.4 


2335. 


830 b 


0.962 


-17.6 


108.7 


2336. 


809 b 


0.134 


71.8 


-180.3 


2337. 


786 b 


0.306 


87.4 


-233.1 


2338. 


770 b 


0.479 


8.8 




2339. 


78 l b 


0.657 


-77.0 


204.0 


2381. 


748 b 


0.044 


22.9 


-68.1 


2382. 


714 b 


0.214 


91.2 


-241.5 


2383. 


712 b 


0.390 


58.9 


-159.7 


3130. 


693 c 


0.865 


-70.1 


185.1 


3131. 


73 l c 


0.048 


26.8 


-70.0 


3132. 


765 c 


0.230 


93.7 


-247.8 


3133. 


800 c 


0.412 


46.7 


-127.4 


3134. 


675 c 


0.566 


-38.9 


97.5 


3135. 


718 c 


0.750 


-93.9 


245.4 



Note. — Instrumental setups: a ESO CAT + CES; b ESO 1.5m + FEROS; c ESO/MPG 2.2m + FEROS. 
RVs are quoted in the zero velocity reference frame. For future conversion, we adopted the systemic velocities 
7T= -21.5kms- 1 and ^ = -23.1kms~ 1 (see 



Table 2: Journal of the XMM-Newton observations. 



JD He i Effective duration (ks) Count rates (10 3 cnt s 1 ) 



-2 450 000 




MOSl 


MOS2 


pn 


MOSl 


MOS2 


pn 


2158.214 


0.700 


33.1 


33.2 


30.6 


10.1 ±0.7 


11.4 ±0.7 


36.0 ± 1.4 


2158.931 


0.826 


19.8 


19.8 


16.5 


10.7 ± 1.0 


10.4 ±0.9 


40.2 ±2.1 


2159.796 


0.978 


33.7 


33.9 


30.1 


11.0 ±0.7 


9.8 ±0.7 


36.7 ± 1.5 


2160.925 


0.177 


26.0 


24.3 


11.7 


12.0 ±0.9 


11.3 ±0.9 


46.0 ±2.6 


2161.774 


0.326 


30.9 


31.0 


28.4 


11.9 ±0.8 


12.5 ±0.8 


38.6 ±1.5 


2162.726 


0.494 


32.9 


32.8 


30.3 


9.8 ±0.7 


10.6 ±0.7 


36.0 ±1.4 



Table 3: Orbital solutions deduced from different RV data sets (i.e. associated to different 
lines). 



Lines 


P 




1 ic 


K, 


K 2 


7i 


72 


r.m.s. 




(d) 




HJD-2 450 000 


(kms -1 ) 


(kms -1 ) 


(kms -1 ) 


(kms -1 ) 


(kms -1 ) 


He I A4471 


5.68151 


4.7 


3199.634 ± 0.004 


94.2 ± 0.5 


245.5 ± 1.2 


-22.6 ± 0.4 


-23.0 ± 0.5 


2.0 


He I A4922 


5.68149 


4.2 


3199.627 ± 0.006 


94.5 ± 0.8 


245.7 ± 2.0 


-20.7 ± 0.7 


-22.1 ± 0.9 


3.2 


HelA50f6 


5.68153 


3.7 


3199.640 ± 0.003 


93.3 ± 0.4 


247.9 ± 1.0 


-22.5 ± 0.3 


-26.4 ± 0.3 


1.9 


He I A5048 


5.68150 


2.1 


3199.639 ± 0.004 


92.4 ± 0.7 


249.7 ± 1.7 


-22.4 ± 0.6 


-25.4 ± 0.7 


3.2 


HelA5876 


5.68148 


3.6 


3199.630 ± 0.002 


93.2 ± 0.3 


246.7 ± 0.7 


-20.5 ± 0.3 


-20.8 ± 0.3 


1.3 


HelA7065 


5.68147 


3.4 


3199.641 ± 0.003 


94.6 ± 0.4 


242.4 ± 1.0 


-19.6 ± 0.5 


-20.9 ± 0.4 


1.8 


He II A4686 


5.68154 




3199.640 ± 0.002 


94.6 ± 0.5 




-16.5 ± 0.2 




1.5 


IIIA5592 


5.68150 




3199.630 ± 0.002 


94.7 ± 0.3 




-24.4 ± 0.2 




1.1 


C IIIA5696 


5.68155 




3199.643 ± 0.002 


93.8 ± 0.8 




-34.3 ± 0.2 




3.1 



a Ti c is the time of inferior conjunction, the primary being in front of the secondary. 
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Table 4: Orbital and physical parameters of CPD — 41°7733. 

Prim. He I lines Lit. 

P (d) 5.681514 5.681504 5.681534 

±3.9110~ 4 ± 3.18 10" 4 ±4.710~ 5 

Syj ' s x ■ ■ ■ 3.3 • • ■ 

Mi/M 2 ■ • • 2.640 ± 0.012 

T ic (HJD 3199.632 3199.635 3199.640 

-2 450 000) ± 0.002 ± 0.003 ±0.002 

K x (kms- 1 ) 93.6 ± 0.3 93.7 ± 0.3 91.9 ± 0.2 

K 2 (kms -1 ) ■•• 247.4 ± 0.9 

v ,i (kms" 1 ) -0.3± 0.3 0.0 ± 0.2 0.0 ± 0.1 

w ,2 (kms" 1 ) ■•• -0.3 ± 0.3 

a lS ini(R ) 10.51 ± 0.06 10.52 ± 0.06 10.31± 0.05 

a 2 sini(R ) ■•• 27.76 ± 0.10 

M lS in 3 i(M ) ■•• 16.94 ± 0.15 

M 2 sin 3 i(M ) ■•• 6.42 ± 0.05 

/ mass (M ) 0.483 ± 0.004 ■ • • 0.457 ± 0.003 

r.m.s. (kms -1 ) 1.5 1.9 5.7 



Note. — Note that the present fits are done in the zero velocity frame. Hence the systemic velocities in 
this frame (noted vg.i) are expected to be zero (see H3.2l and M3.3I for values in the heliocentric frame). 
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Table 5: Average equivalent widths (W\) of selected CPD — 41°7733 spectral lines 



Line 


W A) i 


w Xj2 




(mA) 


(mA) 


He IA4121 


161: 


21: 


He I A4144 


188 ± 15 


61 ± 12 


He I A4388 


273 ± 21 


83 ±25 


He I A4471 


635 ± 35 


100: 


HelA5016 


255 ±5 


20 ±4 


He I A5876 


945 ± 11 


58 ±2 


He II A4542 


382 ± 28 




He II A4686 


590 ± 20 




C IIA4267 


39 ±3 


11 ± 2 


C IIIA5696 


104 ±8 




N II A4630 


105 ±5 


7: 


II A4070 


blend 


9: 


IIIA5592 


251 ±7 




Si IIIA4553 


71 ±4 


7 ±2: 


Si IIIA4575 


29 ±5: 


7: 


Si IVA4089 


354 ± 16 




Mg II A4481 


95 ±9 


20 ±3 



Note. - 



The W\ are referred to the total flux of the binary system. 



Table 6: EPIC spectra best fit parameters and observed X-ray fluxes. 



He i kTi norm x 7V H , 2 a kT 2 norm 2 xl (d.o.f.) / x /x,s /x,m /x,h 



(keV) (l(r 5 cm- 5 ) (10 22 cm- 2 ) (keV) (lO^cm- 5 ) (l(r 14 ergcm- 2 s" 1 ) 



0.700 


o.268;l 


10.2& 6 


U - O4 0.36 


o.73g : g 


5.8l:| 


1.06 (87) 


7.3 


4.6 


2.4 


0.2 


0.826 


o.i98:?I 


12.4 21 4 4 


26 o n 


o.508:|? 


O c:19.8 
y -°5.8 


0.82 (56) 


8.0 


5.3 


2.7 


0.1 


0.978 


0.28°;| 


9.6& 


o.58°|! 


74 - 83 

u - ' ^0.63 


d - u 3.6 


0.80 (94) 


7.5 


4.8 


2.5 


0.2 


0.177 


o 2 c ; - 30 




0.80 i| 


67 2 - 58 

u -° ' 0.50 


12.3 23 ! 8 


1.23 (54) 


8.2 


4.7 


3.2 


0.3 


0.326 


26 31 


10.4™ 


47 - 91 

u -^' 0.16 


63 112 

U.DO 52 


7 i 10.4 
' • ± 2.2 


1.14 (89) 


7.9 


5.0 


2.7 


0.2 


0.494 


0.28^4 


8.6? : § 


41 °- 74 


71 °- 87 

u - ' i 0.62 


4.8§ : f 


1.25 (87) 


7.0 


4.5 


2.4 


0.2 



Merged 0.28g;l§ 9.9^° 3 5 0.54g;g 0.80j:?l 5.1^ 1.37 (315) 7.7 4.8 2.6 0.2 



a As indicated by its subscript, the local equivalent Hydrogen column 7Vh,2 only applies to the second tem- 
perature component kT 2 of the model. 
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Table 7: Absorption corrected fluxes and X-ray luminosity of CPD — 41°7733. 



0He I 


-fun. 

Jx 


fun. 

Jx,s 


fun. 
JX,M 


fun. 
JX,B 


log L x 




(io- 


-14 —2 - 

14 erg cm z s 




(ergs" 1 ) 


2-T models 


0.700 


21.3 


17.6 


3.5 


0.2 


31.83 


0.826 


25.7 


21.6 


3.9 


0.1 


31.92 


0.978 


21.0 


17.1 


3.6 


0.2 


31.83 


0.177 


23.0 


18.3 


4.4 


0.3 


31.87 


0.326 


22.8 


18.8 


3.9 


0.2 


31.86 


0.494 


19.6 


16.0 


3.4 


0.2 


31.80 


Merged 


21.5 


17.5 


3.8 


0.2 


31.84 


3-T model 


Merged 


22.4 


18.0 


3.8 


0.6 


31.86 


2-T + PL model 


Merged 


22.4 


18.0 


3.8 


0.6 


31.86 
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Fig. 1. — Combined epic pn image in the range 0.5-10.0 keV. Adopted source and pn back- 
ground (bkg) extraction regions are indicated. North is up, East to the left. 
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A (A) 

Fig. 2.— Blue spectrum of CPD -41°7733 as obtained on HJD = 2 453 135.718 (0 = 0.750). 
The identifications of the main lines have been indicated. The He I lines display a clear SB2 
signature, with the blue-shifted primary and the red-shifted secondary components. The 
identification ticks refer to the rest wavelengths. 
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Fig. 3. — He I A4922 (left) and He I A5876 (right) lines at selected phases. The spectra were 
shifted along the vertical axis according to the value of their phase (right-hand scale). The 
secondary spectral signature is clearly identified in all but the blended spectra. 
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0.5 1 



Fig. 4. — CPD — 41°7733 RV curves corresponding to the circular He I solution of Table 
HI The RV measurements listed in Table [T] are also displayed. Filled and open symbols 
correspond to the primary and secondary components respectively. The squares give the 
CES data while the circles indicate the FEROS measurements. 
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i r 




i 



5.6 5.65 5.7 5.75 5.8 

P (d) 



Fig. 5. — Root-mean square (r.m.s.) residuals of the best circular orbital solutions as a 
function of the adopted period. The adjust ments we re performed using the different RV sets 
presented in Fig. El excluding the point of IPHYB9G at RV = — 109kms _1 . 
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Fig. 6. — Combined RV curve computed using literature data together with the measure- 
ments of Table [TJ Diff erent symbol s refer to different data sets. Open symbols are for 
previous obser vations : Struve 19441. a sterisk; PHYB90 . upward triangles; HCB74 . down- 
ward triangles; ILM83I . diamonds; ISLOll . square. Filled symb ols indica te new RV points from 
this work: CES, squares; FEROS, circles. Note that the IpHYB9oI point at RV = -109 
kms -1 , near « 0.5, was excluded for the period search and orbital solution determination. 
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4. 



4. 



4 



eff 



Fig. 7. — Location of the CPD — 41°7733 primary (filled symbols) and secondary (open 
symbols) components in the H-R diagram. The triangles (resp. squares) indic ate an adopted 



giant (resp. main sequence) luminosity class. The evolutionary tracks from ISchaller et al. 



(119921 ) are shown (plain lines) as well as the isochrones (dotted lines) computed for ages 
ranging from 2 to 10 Myr with a step of 2 Myr. The boxes present the ranges of parameter 
values obtained assuming different luminosity classes for the components. 
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Fig. 8. — Merged EPIC MOSl (black), MOS2 (red) and pn (green) spectra of CPD -41°7733 
simultaneously fitted with 2-T (upper panel) and 3-T (lower panel) mekal models. The 
bottom window of each panel shows the contributions of individual bins to the \ 2 of the fit. 
The contributions are carried over with the sign of the deviation (in the sense data minus 
model). The figures appear in color in the electronic version of the paper. 



